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In this study, two structural models (denoted as M’’°ST and M?°SH) of the main proteinase (M?’) from 
the novel coronavirus associated with severe acute respiratory syndrome (SARS-CoV) were constructed 
based on the crystallographic structures of M”® from transmissible gastroenteritis coronavirus (TGEV) 
(MPT) and human coronavirus HcoV-229E (M?’H), respectively. Various 200 ps molecular dynamics simu- 
lations were subsequently performed to investigate the dynamics behaviors of several structural features. 
Both M?°ST and MSH exhibit similar folds as their respective template proteins. These structural models 
reveal three distinct functional domains as well as an intervening loop connecting domains II and III as found 
in both template proteins. In addition, domain III of these structures exhibits the least secondary structural 
conservation. A catalytic cleft containing the substrate binding subsites S1 and the $2 between domains I and 
II are also observed in these structural models. Although these structures share many common features, the 
most significant difference occurs at the $2 subsite, where the amino acid residues lining up this subsite are 
least conserved. It may be a critical challenge for designing anti-SARS drugs by simply screening the known 


889 


database of proteinase inhibitors. 
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INTRODUCTION 


An outbreak of atypical pneumonia, designated as se- 
vere acute respiratory syndrome (SARS), was first reported 
in Guangdong Province of China in late 2002, and rapidly 
spread to several countries.'” Infection by SARS is usually 
characterized by high fever, malaise, rigor, headache, non- 
productive cough and may progress to generalized, intersti- 
tial infiltrates in the lung.’ Attempts to identify the etiology 
of the SARS outbreak were not successful until March 2003, 
when laboratories in the United States, Canada, Germany, 
and Hong Kong isolated a novel coronavirus (SARS-CoV) 
from SARS patients. The sequence of the complete genome 
of SARS-CoV was further determined and characterized with 
two different isolates.*° Phylogenetic analyses and sequence 
comparisons reveal that SARS-CoV is not closely related to 
any of the three groups of coronaviruses, including two hu- 
man coronaviruses, HcoV-229E (group I) and HcoV-OC43 
(group II), which are responsible for about 30% of mild upper 


respiratory tract illnesses,°* in particular, the common cold.’ 

Coronaviruses belong to a diverse group of positive- 
stranded RNA viruses featuring the largest viral RNA genomes 
known to date (27-31 kb). They share a similar genome orga- 
nization and common transcriptional and translational pro- 
cesses as Arteriviridae.'°'' The human coronavirus HcoV- 
229E replicase gene encodes two overlapping polyproteins, 
ppla (replicase la, ~450 kDa) and pplab (replicase Lab, 
~750 kDa),!” that mediate all the functions required for viral 
replication and transcription.'* The functional polypeptides 
are released from the polyproteins by extensive proteolytic 
processing, which is primarily achieved by the 33.1-kDa 
HCoV-229E main proteinase (M””’).'* M”” is commonly also 
called 3C-like proteinase (3CL"”) to indicate a similarity of 
its cleavage site specificity to that observed for picornavirus 
3C proteinase (3C”) and the identification of a Cys residue 
as the principle nucleophile in the context of a predicted 
two-B-barrel fold.'*'® M’° from HcoV-229E (M?H) has 
been biosynthesized in Escherichia coli and the enzyme 
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properties, inhibitor profile, and substrate specificity of the 
purified protein have been well characterized.'*'” 

Several studies have revealed significant differences in 
both the active sites and domain structures of M?® from coro- 
navirus and picornavirus.'’*! It is noteworthy that corona- 
virus M"® possesses a large C-terminal domain of ~110 
amino acid residues (domain III) that is not found in other 
RNA virus 3CL?”. Deletion of this domain results in dra- 
matic losses of proteolytic activity, suggesting that the 
C-terminal domain III of M’” contributes to proteolytic ac- 
tivity through undefined mechanisms. Previous experimental 
data have shown that the differential cleavage kinetics of 
sites within ppla/pplab are a conserved feature of corona- 


pro 


virus M’® and that similar processing kinetics for the repli- 
case polyproteins of all coronaviruses can be predicted.” 
Furthermore, the cleavage pattern appears to be conserved in 
MP? from SARS-CoV (M?°S) and from other coronavir 
uses,” as deduced from the genome sequence.” The func- 
tional importance of M”” in the viral life cycle has made this 
proteinase an attractive target for the development of drugs 
directed against SARS and other coronavirus infections. 
Moreover, molecular modeling has suggested that available 
rhinovirus 3C?" inhibitors such as compound AG7088 may 
be modified to be tested for SARS therapy.” Therefore, 
screening the known proteinase inhibitor libraries may be an 
appreciated shortcut to discover anti-SARS drugs.” 
Recently, crystal structures of MH” and M”” from 
porcine coronavirus (transmissible gastroenteritis virus, 
TGEV) (M?°T) complexed with its inhibitor” have been de- 
termined. In addition, homology models of M?S based on 
the crystal structures of MH” and M?°T” have been also 
reported. Comparison of these structures reveals a remark- 
able degree of conservation of the substrate binding sites, 
which is further supported by the cleavage of the substrate for 
the M?°T with the recombinant M"’’S.”* In addition, M?"°S 
shows 40 and 44% sequence identity to MP’°H and M?°T, re- 
spectively.”* Although the results from the deduced genome 
sequence of SARS-CoV have indicated that it belongs to a 
new group of coronaviruses,” the significantly high se- 
quence identity of M’"S to bovine coronavirus (BCoV) M?® 
(49%) and mouse hepatitis virus (MHV) M? (50%) from 
group II coronaviruses has allowed Anand et al.”? to recog- 
nize it as an outlier among group II coronaviruses. 
Molecular dynamics (MD) simulations in the atomic 
level have been intensively applied to gain insight into the 
structure-function relationships of proteins. Previously, sev- 
eral MD simulations and molecular docking experiments 
have been successfully conducted towards various target pro- 
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teins in our group.”’* In this paper, two structural models of 


MS (denoted as M’'°SH and M?'°ST) were constructed 
based on the crystallographic structures of M?’H?* and 
M?°T”°, respectively, by the comparative approach. In addi- 
tion, MD simulations were conducted to investigate the dy- 
namics behaviors of these structures. Beyond the continued 


pro 


characterization of M”” from various coronaviruses, the 
amino acid sequence alignment and structural homology 
analyses of M?'’S presented in this study provide particularly 
attractive targets for further structure-based studies, such as 
folding/unfolding mechanism and molecular docking, which 


are currently being carried out in our group. 


METHODS 


Model proteins 

Structural homology to construct the structural models 
of MPS (MST and M?°SH) was based on the monomer of 
the three-dimensional (3D) structure of M?°T, refined to 1.96 
A resolution”® (Fig. 1(A)), and that of MH, solved at 2.54 A 
resolution” (Fig. 1(B)), obtained from the protein data bank 
(PDB; accession numbers llvo and 1p9u, respectively). The 
inhibitor, a substrate analog hexapeptidyl chloromethyl ke- 
tone, was removed from the crystallographic structure of 
M?"’T before being used as a template. Unfavorable nonphys- 
ical contacts in these structures were then eliminated using 
the Biopolymer module of the Insight II program (Accelyrs, 
San Diego, CA, USA) with the force field Discover CVFF 
(consistent valence force field)**** in the SGI 0200 worksta- 
tion with a 64-bit HIPS RISC R12000 2 x 270 MHz CPU and 
PMC-Sierra RM7000A 350 MHz processor (Silicon Graph- 
ics, Inc., Mountain View, CA, USA), followed by 10,000 en- 
ergy minimization calculations using steepest descent method, 
to yield the model proteins for further structure building. 


Structural homology 

Homology utilizes structure and sequence similarities 
for predicting unknown protein structures. The Homology 
module in Insight II allows us to build the 3D models of the 
target protein (i.e., M”°S) using both its amino acid sequence 
and the structures of known, related model proteins (i.e., 
MPH and MT). The Homology program provides simulta- 
neous optimization of both structure and sequence homo- 
logies for multiple proteins in a 3D graphics environment, 
based on a method developed by Greer.*’ Smith-Waterman 
pairwise amino acid sequence alignments were performed 
based on the conserved active site and substrate binding 
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subsites among M?” from various coronaviruses to find the the target protein were generated from alignments of the tar- 
location of the active site and substrate binding subsite of — get protein to the model proteins. The atomic coordinates 
MPS. The consensus structural conserved regions (SCRs) of _ were then transferred from the model proteins to the target 


Fig. 1. The x-ray crystallographic structure of (A) M?'°T and (B) M?'°H and the structural model of (C) M?*°ST and (D) 
M?"°SH. These structures are visualized by the Insight II program. The N- and C-termini are indicated. o-Helices are 
shown in red cylinders, while B-strands are illustrated in arrows pointing from N- to C-terminus. The polypeptide 
backbones belonging to the turn and random coil regions are shown in blue and green, respectively. The general acid- 
base catalyst His residue and the nucleophilic Cys residue are labeled. The locations of the putative substrate binding 
subsites S1 and S2 are indicated. 


892 J. Chin. Chem. Soc., Vol. 51, No. 5A, 2004 


protein in each SCR using the Mutation Matrix module of the 
Insight II program. Automatic loop building was performed 
either by database searching” or generation through random 
conformational search.*’ The coordinates at the N- and C- 
termini of these loops were then automatically assigned. Side 
chains of the target protein were automatically replaced, pre- 
serving the conformations of the model proteins. The side 
chain conformations were optimized either manually or auto- 
matically using a rotamer library.” Similar secondary struc- 
ture motifs were identified by database searching and pre- 
dicted by DSSP.*° The bond lengths and torsion angles in the 
SCRs and loop regions were repaired and relaxed using 
Homology/Refine/SpliceRepair and Homology/Refine/Re- 
lax, respectively. The newly built structures of the target pro- 
tein were substantially refined to avoid van der Waals radius 
overlapping, unfavorable atomic distances, and undesirable 
torsion angles using molecular mechanics and dynamics fea- 
tures in the Discover module. 


Molecular dynamics simulations 

The crystallographic structures of M?°H and M?°T and 
the structural models of M?’°SH and M?°ST were subjected to 
energy minimization calculations by steepest descent method 
with 3,000 iterations followed by Newton-Raphson method 
with 5,000 iterations to be used as the initial energy-mini- 
mized structures for further structural comparison. Each en- 
ergy-minimized structure was subsequently placed in the 
center of a lattice with the size of 50 x 60 x 85 A® full of 
6,222, 5,866, 5,836, and 5,776 water molecules for the sys- 
tem of M?°H, MT, M?°SH, and M?'°ST, respectively. 
These systems composed of the target protein and water mol- 
ecules were then equilibrated by performing 20,000 steepest 
descent minimization and 10 ps dynamics calculations. The 
explicit image periodic boundary condition (PBC) was used 
for solvent equilibrium. At the end of explicit image equilib- 
rium, Discover will re-image a molecule whose center of 
mass has moved out of the lattice in order to maintain the in- 
tegrity of the lattice with a relatively constant density. A cut- 
off of 14 A was used to calculate long-range electrostatic in- 
teractions. Finally, 200 ps MD simulation was carried out for 
each system using the Discover module of Insight II. The 
temperature and pressure were maintained for each MD sim- 
ulation by weak coupling the system to a heat bath at 300 K 
and an external pressure bath at one atmosphere with a cou- 
pling constant of 0.5 ps, according to the method described 
by Berendsen et al.*! A cut-off radius of 10 A for the non- 
bonded interactions was applied to each MD simulation. The 
time-step of the MD simulations was | fs. The trajectories 
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and coordinates of these structures were recorded every 2 ps 
for further analysis. 


Structural analyses 
Although some complicated algorithms have been pro- 
posed to measure the structural similarity between pro 


s 42,43 
teins,” 


the root-mean-square deviation (RMSD) remains 
the simplest one for closely related proteins.** For each MD 
simulation, the RMSDs of the trajectories recorded every 2 ps 
interval were calculated for the backbone C, atom of the en- 
tire protein, domains J, II, and II, and the substrate binding 
subsites $1 and S2 during the course of 200 ps MD simula- 
tions with reference to the respective starting structure ac- 
cording to Koehi.*” The RMSDs were calculated after opti- 
mal superimposition of the coordinates to remove transla- 
tional and rotational motion.*° Secondary structures were 
predicted based on DSSP,“° in which pattern recognition of 
the hydrogen bond was correlated to the geometrical features. 
The default hydrogen bonding energy criterion of -0.5 kcal/ 
mol was used. Accessible surface areas (ASAs) of the sub- 
strate binding subsites $1 and S2 and the linear distance be- 
tween the sulfur atom of the nucleophilic Cys residue and the 
N® of the general acid-base catalyst His residue for each 
structure were also recorded as a function of MD simulation 
time. 


RESULTS AND DISCUSSION 


Amino acid sequence alignment 

Point, insertion, or deletion mutations that would result 
in a critical loss of biological functions are less favored by 
evolution and consequently, functionally and structurally rel- 
evant domains tend to be highly conserved across a corre- 
sponding protein family. Such conservation can be detected 
as a pattern of conserved residues that would be unlikely to 
have occurred by chance. Therefore, an optimal amino acid 
sequence alignment based on the conserved residues is essen- 
tial to the success of structural homology. The results of 
amino acid sequence alignment of M"°S to M?°T and M?°H 
are given in Fig. 2. There are 301, 300, and 306 residues in 
MPT, MPH, and MS, respectively. The residue corre- 
sponding to Ala46 in domain I of M?'°S and those corre- 
sponding to Asp248, e249, and Gln273 in domain III of 
MPS are missing in both M?°T and MPH. In addition, there 
are one and two extra residues at the C-terminus of M?°S 
compared to M?°T and MP"°H, respectively. There are 33, 33, 
and 29 totally conserved residues in domains I, II, and III 
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among these MP", respectively, indicating that domain III ex- 
hibits higher sequence variation among these three domains. 
It has been predicted earlier that the active site of the corona- 
virus M”” is similar to those of picornavirus 3C””.'° Both the 
general acid-base catalyst and the nucleophile residues of 
these three proteins are totally conserved, where the general 
acid-base catalyst His41 is located in a highly conserved sig- 
nature sequence (LNGLWLXDXVXCPRHVJ) of domain I 
and the nucleophilic Cys144 for MP°T and M?°H or Cys145 
for M?'°S is located in the highly conserved signature se- 
quence (TIXGSFXXGXCGSXG) of domain II (i.e., Xs indi- 
cate the nonconserved residues). 

Table | lists the percentages of amino acid identity 
among these proteins. M?°T and M”’°H show the highest total 
amino acid identity (60.80%), whereas M?°H and MPS ex- 
hibit the lowest total amino acid identity (40.19%). In addi- 
tion, domain II has the highest amino acid identity, whereas 
domain III shows the lowest amino acid identity among these 
three proteins. M’'°S shows slightly higher amino acid iden- 
tity to MPT than M?"H, indicating that the structure of MPS 
may be more similar to M?°T than M?°H. The above results 
are in good agreement with the previous finding that MS 
shows 40 and 44% sequence identity to M”°H and M?"T, re- 
spectively.”* Although the significantly high sequence iden- 
tity of MPS to BCoV MP” (49%) and MHV M?”” (50%) from 
group II coronaviruses has allowed Anand et al.”* to recog- 
nize it as an outlier among group II coronaviruses, the low se- 
quence identities between M?°S and M?°T and between 
MPS and M?°H from the present study strongly support the 
results of Marra et al.,* in which SARS-CoV was classified as 


a new group of coronavirus from the analysis of the deduced 
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Table 1. The amino acid sequence identities among M° H, 
MPT, and MPS 


Identity (%) 
Total Domain I Domain II Domain III 
MPH and MP°T 60.80 63.44 65.06 55.45 
MPH and MP°S ~—s 40.19 41.94 45.78 35.64 
MPT and MP?S) 43.85 44.09 49.40 39.22 


genome sequence. 


The Structural Models of M?°ST and M?'SH 

Recently, two 3D models of M?'°S have been con- 
structed based on the crystallographic structure of MH”? 
and that of M’°T,”* using the homology modeling technique. 
Virtual screening was further performed employing molecu- 
lar docking towards both constructed models to identify pos- 
sible 3CL?" inhibitors from small molecular databases for 
SARS therapy. The level of similarity between M?°S and 
M?°T as well as M?°H allowed us to construct two structural 
models for M?°S (denoted as M?'°ST and M?'°SH) by the 
comparative approach and the results are illustrated in Figs. 
1(C) and 1(D). There are three 1- and 2-residue insertions in 
MS, relative to both structural templates (Fig. 2); as to be 
expected, these are all located in loops and do not present a 
problem in model building. Both M?°ST and M?°SH exhibit 
three distinct domains, indicating that they adopt similar 
folds as MP°T and M”°H, respectively. However, the second- 
ary structures of both M’°ST and M?SH predicted accord- 
ing to DSSP*° were less conserved compared to those of 
MPT (Fig. 1(A)) and M?°H (Fig. 1(B)), particularly in do- 


Fig. 2. Amino acid sequence alignment of M?°T, MP? H, and MP"°S. Secondary structures as defined in the crystallographic 
structure of M?'°T are shown on top. The start and end amino acid residues are numbered in the brackets on the left and 


right of each sequence, respectively. Residues totally conserved in all sequences are indicated in red letters with green 
background. Residues conserved in M?°T and M?'°H but different from those in M?"°S are represented in black letters 
with yellow background. Residues where variations occur are given in blue or brown letters with grey background. 
The amino acid residues missing in both M?°T and M?'°H are shown as dashed lines. 
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main III. The result is consistent with that of amino acid se- 
quence alignment, showing that domain III exhibits the least 
sequence identity compared to domains I and II among these 
proteins. It further implies that some of the main-chain or 
side-chain hydrogen bonds in the constructed homology 
models may be destroyed in order to maintain folds similar to 
the model proteins. 

The putative substrate binding subsites $1 and S2 of 
M?ST and MSH are located in a cleft between domains I 
and II, which are nearly identical to those of M?'°T and M?°H 
(Fig. 1). It indicates that M?°S may follow the similar sub- 
strate binding mechanisms of M?°T and M?°H, allowing us 
to design anti-SARS drugs by screening the known pro- 
teinase inhibitors. Instead of separating domains I and II with 
the catalytic cleft, a long loop (residues 184 to 199 in both 
M?°T and M?'°H and residues 185 to 200 in M?°S) loosely 
connect domains II and III in all structures. Domain III, a 
globular cluster of 5,5, 4, and 2 helices for M?'°T, M?'°H, 
MP?PST, and M’°SH, respectively (Fig. 1), has been impli- 
cated in the proteolytic activity of M"°.'’ Indeed, there 
should be only 3 and 1 helices in domain HI of M?°ST and 
M°SH, respectively, both with helix AIL broken into two 
parts. Comparing the two crystallographic structures, MP°T 
and M?'°H, and the two homology models, M”°ST and 
MSH, we found that domain I of M?"°S is more similar to 
that of M?'°H, while domains II and HI of M’’’S are more sim- 
ilar to those of M?°T. The low sequence identity and second- 
ary structure similarity in domain III among these proteins 
presented in the present study, as well as the previous find- 
ings showing that the characterization of recombinant pro- 
teins, in which 33, 28, and 34 C-terminal amino acid residues 
of MP? from IBV, MHV, and HCoV, respectively, were de- 
leted resulted in dramatic losses of proteolytic activity, sug- 
gest that domain II may play a minor role in proteolytic ac- 
tivity through an undefined mechanism. 

The analysis of RMSD (Table 2) shows that the struc- 
ture of M’'°H is very similar to that of M?°T,?° with the 
RMSD between these two structures being 2.01 A for all 300 
C. positions of the molecule. M?'°H used in this work lacks 
two amino acid residues from the C-terminus. Nevertheless, 
it has the same enzymatic properties as full-length M’°H but 
yields much superior crystals.” In the structure of full-length 
M?"°H, residues 301 and 302 are disordered and not seen in 
the electron density.”’ With both HCoV 229E and TGEV be- 
ing group I coronavirus,”’ their M”° share 60.8% sequence 
identity (Table 1). As shown in Table 2, the RMSDs of the 
constructed models, MSH and M”°ST, are 4.84 and 3.94 A, 
compared to their corresponding templates, M?’°H and M?°T, 
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Table 2. The RMSDs between the template proteins, M?'°H and 
MPT, and the structural models, M?°SH and M?"°ST 


RMSD (A) 
MP°H MP°T M?°SH 
MPT 2.01 Z = 
MP?SH 4.51 3.94 < 
MP°ST 4.84 4.37 5.78 


respectively, while the RMSD between M”SH and M?"’ST is 
5.78 A. It indicates that the structure of M’’S is more similar 
to that of M?°T than that of M?'° H. It further implies that the 
choice of a more closely related template protein yields a 
more accurate structural model of MPS. 


Molecular dynamics simulations 

The six monomers of M?°T” and M’°H”* presented the 
asymmetric unit are arranged in three dimmers. Each mono- 
mer is folded into three distinct domains, the first two of 
which are antiparallel §-barrels reminiscent of those found in 
serine proteinases of the chymotrypsin family. Residues 
8-100 form domain I, and residues 101-183 make up domain 
II (Fig. 2). The connection to domain III is formed by a long 
loop comprising residues 184-199. Domain III, composed of 
residues 200-302 (Fig. 2), contains a novel arrangement of 5 
a-helices (Fig. 1(A) and (B)). A deep cleft between domains I 
and II, lined up by hydrophobic residues, constitutes the sub- 
strate binding subsites. The catalytic site is situated at the 
center of the cleft. In order to investigate the dynamics behav- 
iors of MPT, MPH, MP'°ST, and M?'°SH, various MD simu- 
lations of these proteins in explicit water were conducted at 
300 K. The overall structural changes were evaluated by plot- 
ting the RMSDs of these proteins relative to the original posi- 
tions in the corresponding starting structures as a function of 
running time as shown in Fig. 3. During the MD time course, 
these structures remained considerably stable, with the 
RMSDs remaining within 3 A. It is obvious that domain III 
exhibited more structural variations than the other two do- 
mains in all cases. The substrate binding subsite S1 of M?°ST 
was found to maintain its structural integrity during the entire 
MD time course, whereas the substrate binding subsite S2 of 
M?’ST exhibited more structural variations. The higher 
structural variation of S2 makes it flexible enough to accom- 
modate a bulky hydrophobic residue from the substrate. 

The secondary structure propensity of these proteins 
was predicted according to DSSP*° during the entire MD 
courses and the results are shown in Fig. 4. As expected, both 
M?T and M?°H exhibited higher secondary structure stabil- 
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ity than M’'°ST and M”°SH. The interior of the B-barrel of 
domain I consists entirely of hydrophobic residues. A short 
a-helix (Helix AI; Tyr53-Ser58) closes the barrel like a 
lid.??”° However, this short a-helix in both MT and MH 
is unstable during the 200 ps MD simulations (Fig. 4(A) and 
(B)). Furthermore, this short a-helix is missing in M?°ST and 
MPSH (Fig. 4(C) and (D)). Domain I is smaller than domain 
I and also smaller than the homologous domain II of shy- 
motrypsin and hepatitis A virus (HAV) 3C”".**°° Several 
secondary structure elements of HAV 3C”” are missing in 
both M?°T and M?°H. Domain III is composed of 5, mostly 
antiparallel, a-helices and the loops connecting them. Inter- 
helical contacts are mediated by hydrophobic side chains.”° 
Database searches’”! did not reveal other proteins or protein 
domains with the same topology as domain III. The present 
homology models showed that some of the secondary struc- 
tures of M?°T and M?°H were missing in domain III. It is 
possible that in order to maintain folds similar to the template 
proteins, some of the main-chain and side-chain hydrogen 
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bonding patterns of both MST and M?°SH were missing, 
resulting in the loss of the secondary structure content. Our 
results again suggest that domain III of these M”” play a role 
in proteolytic activity through an undefined mechanism re- 


gardless of its structural integrity. 


Active site 

The active site of the coronavirus M”” is similar to 
those of the piconavirus 3C””, as had been predicted previ- 
ously.'° The mutual arrangement of the nucleophilic Cys144 
and the general acid-base catalyst His41 of M?'°T is identical 
to that of the HAV 3C”® Cys172 and His44 and the Ser195 
and His57 residues of chymotrypsin.” The distance between 
the sulfur atom of Cys144 and the N® of His41 in M?”T is 
4.05 A,”* longer than the corresponding Cys-His distances in 
HAV 3C?? (3.92 A), poliovirus (PV) 3C”® (3.4 A), and 
papain (3.65 A).** From a dynamics point of view (Fig. 5), the 
Csy144-His41 distance of M?°H fluctuated more rapidly 
than that of M’’°T. In addition, the Cys145-His41 distances of 


Time (ps) 


Fig. 3. The RMSDs of the backbone C, for (A) the whole protein, (B) domain I, (C) domain II, (D) domain III, (E) substrate 
binding subsite S1, and (F) substrate binding subsite S2 of MPT, M?*°H, MP*°ST, and MP?'°SH with reference to their 
respective starting structure during the 200 ps MD simulations at 300 K. 
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MP?SH fluctuated more rapidly than that of MST beyond 
150 ps. These results indicate that both M?°T and M?°ST 
may exhibit more stable active site configurations than those 
of MPS and MP°SH. 

It is generally assumed that the native state of the active 
site of papain-like cystein proteinases is a thiolate-imidazo- 
lium ion pair formed by Cys and His residues.” In protein- 
ases of the papain family, an Asn residue is the third member 
of the catalytic triad. Chymotrypsin and other members of 
this serine proteinase family have a catalytic triad consisting 
of Ser195-His57-Asp102. In HAV 3C?", Asp84 is present at 
the required position.””° PV 3C”", human rhinovirus (HRV) 
3C”, and HRV 2A?” have a Glu or Asp residue in the proper 
orientation to accept a hydrogen bond from the active site His 
residue.”°°*” Both M?°T and M?°H have Val84 in the corre- 
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sponding position, with its side chain pointing away from the 
catalytic site. The corresponding residue in M?'°S is Cys85 
(Fig. 2). In both M?°T and MP°H, the polypeptide segment 
184-199, which connects domains II and III and is probably 
involved in substrate binding, is held in the proper position 
during catalysis. The corresponding segment was also found 
in MS, although its amino acid sequence is not conserved 
compared to those of M?°T and MPH (Fig. 2). A direct in- 
volvement of His163 or Asp186 of M?°T, Gln163 or Asp186 
of M?° H, and His164 or Asp187 of M?S in catalysis, makes 
them a clear case of viral cystein proteinase employing only a 
catalytic dyad.”° 


Substrate-binding subsites 


It has been shown previously that, similarly to 
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tively. 
3cP”,°?5° specific substrate binding by M”” is ensured by 
well-defined $1 and S2 binding pockets.” In addition, it has 
also been shown that the imidazole side chain of a conserved 
His residue, which is located in the center of a hydrophobic 
pocket, interacts with the Pl carboxamide side chain of the 
substrate. This specific interaction is generally considered to 
determine the piconavirus 3C"” specificity for Gln residue at 


MPH or His163 of MS is located at the very bottom of this 
hydrophobic pocket which is formed by the totally conserved 
residues Phe139 of both M?°T and MP°H or Phe140 of M?°S 
and the main-chain atoms of Ile140, Leul64, Glul65, and 
His171 of M?’T, Ile 140, e164, Glul165, His171 of M?'°H, or 
Leul41, Met165, Glul66, and His172 of MS. The totally 
conserved Glu165 of M?°T and M”°H or Glul66 of M?°S 
forms an ion pair with the totally conserved His171 of MP°T 
and M?°H or His172 of M’°S.”° This salt bridge is itself on 
the periphery of these molecules, forming part of the outer 
wall of the S1 subsite. 

Coronavirus MP?" has a strong preference for Leu resi- 
due at the P2 position.”! Similar to $1 subsite, the putative S2 


subsute identified in the structure is also a hydrophobic 
pocket that is suitably positioned and large enough to accom- 
modate a Leu side chain easily. In both M?°T and M?°H, the 
S2 pocket is lined by the side chains of His41, Thr47, [He51, 
Leul64, and Pro188, despite residue Leu164 in M?°T being 
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replaced by Ie. In M?°S, the S2 pocket is lined by the side 
chains of His41, Asp48, Pro52, Met165, and Gln189. It indi- 
cates that the S2 pocket is not as conserved as the S1 pocket, 
probably because the S2 subsite is located at the open mouth 
of the catalytic cleft formed by domains I and II, while the $1 
subsite is situated in the very bottom of this cleft. Thus, the 
structural conservation of the substrate binding subsite S1 is 
necessary to maintain the structural integrity of both domains 
I and II. It is worth mentioning that the main chain of Leul 64 
of M°T (or He164 of M’°H or Met165 of M?"°S) forms part 
of the S1 subsite while its side chain is involved in the S2 
subsite. It indicates that these two subsites are somewhat in- 
fluenced by each other towards substrate binding. 

Fig. 6 shows the ASAs of both S1 and S2 subsites for 
MT, MPH, MP'ST, and M’'°SH as a function of MD simu- 
lation time. Both subsites are flexible enough to accommo- 
date the substrates. In order to gain a clearer look of these 
structures, the snapshots of both S1 and 82 subsites for these 
proteins with the smallest and largest ASAs sampled from the 
200 ps MD simulations are illustrated in Fig. 7. The smallest 
ASA of S1 is 212, 95.4, 150, and 233 A? sampled at 38, 116, 
146, and 2 ps, while the largest ASA of S1 is 360, 158, 286, 
and 361 A? sampled at 88, 30, 28, and 94 ps for M?°T, MP°H, 
MP®ST, and M?°SH, respectively. The smallest ASA of S2 is 
117, 107, 290, and 143 A’ sampled at 2, 118, 176, and 4 ps, 


400 
350 
300 
250 
200 
al Chr ee Take ee ee mae 
aie »” i oi 
50 

500 

450 
400 
350 
300 
250 
200 LV. 
150 Fe mw eels 
100 L_—__ 


0 50 100 150 200 
Time (ps) 


. 
wa,” 
irae 


ASA (A) 


Fig. 6. The ASAs of the substrate binding subsites (A) 
S1 and (B) S2 as a function of MD simulation 
time for MPT, MP'°H, MP°ST, and MP'°SH. 


Liu et al. 


while the largest ASA of S2 is 344, 217, 461, and 296 A’ sam- 
pled at 158, 196, 26, and 158 ps for M?°T, MP°H, MP°ST, 
and M?'°SH, respectively. It is interesting that the sizes and 
conformations of the smallest and the largest S1 pocket of 
M°SH are very similar to those of M?'°T. The variation of 
the size and conformation of S2 subsite for these proteins is 
more significant than the S1 subsite during the MD simula- 
tions, probably because S2 is fully exposed to the solvent and 
is easy to be subjected to structural change. The structural 
variation of these two subsites allows them to accommodate 
the specific recognition residues of the substrates upon bind- 
ing. 

In conclusion, the technique of the comparative ap- 
proach was successfully applied to construct the homology 
models of M?°ST and MSH based on the crystal structures 
of MT and MP°H, respectively, in this study. Molecular dy- 
namics simulations were subsequently conducted to investi- 
gate the dynamics behaviors of the structural elements of 
these structures. Although these structures share many com- 
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Fig. 7. The snapshots of the substrate binding subsites 
S1 and S2 for MPT, MPH, MPPST, and 
M?'°SH with the smallest and the largest ASAs 
during the 200 ps MD simulations. The protein 
residues are illustrated in CPK with the residues 
forming these subsites being shown in red. The 
residues lining up the hydrophobic pockets of 
subsites S1 and S2 are labeled. The value of the 
smallest and the largest ASAs for each protein 
and the time point the structure was sampled are 
given at the upper-right and the lower-right cor- 
ners of each frame, respectively. 
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mon features, the most significant difference occurs at the $2 
subsite, where the amino acid residues lining up this subsite 
are least conserved. It may be a critical challenge for design- 
ing anti-SARS drugs by simply screening the known data- 
base of proteinase inhibitors. 
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